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1. Introduction 

The universe a few microseconds after the Big Bang was filled with a hot and dense phase 
of matter. We believe that quarks and gluon at those temperatures, above 10^^ K, were de- 
confined and existed as a quark gluon plasma (QGP). These ideas can be tested in collisions 
of nuclei at ultra-relativistic energies. At the Relativistic Heavy Ion Collider (RHIC) nuclei 
as heavy as gold are accelerated to an energy of 100 GeV per nucleon. 40 TeV total energy 
is available in the collision of two of these nuclei. A large fraction of this energy is used to 
create new particles (a few thousand of them) and to give them kinetic energies. The en- 
ergy in the collision is spread over a certain volume with peak energy densities of ^ 10-30 
GeV/fm'^, much larger than in ordinary nuclear matter [1-3]. This process happens over a 
very short time (a few fm/c) but we have strong indications from observables like elliptic 
flow V2 that despite the short life times thermalization of the matter is reached. The peak 
energy densities are safely above the predicted critical energy densities w 1 GeV/fm'^ for 
a phase transition from hadrons to a quark gluon plasma [4]. 

In the first decade of RHIC a large amount of experimental results has been gathered. 
Some results have confirmed expectations, like jet quenching. Others have turned out to 
be more diffult to decipher, e.g. signals from quarkonia that promised to be simple plasma 
thermometers [5], yet others have been great surprises, like the indications for a perfect 
fluid created at RHIC. We do have firm evidence that we create equilibrated quark-gluon 
matter in the experiments at RHIC. [1,2,6] 

In order to describe nuclear collisions at high energies and to filter out the properties of 
quark gluon plasma from the stream of experimental data we can slice the data in a simple 
way. Thermalization of particles does not extend beyond 2 GeV/c in transverse momentum 
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Pt around midrapidities, i.e. around 90° scattering angle from the beam axis. However, 
this includes 99% of all particles produced, the bulk of the collision. The tail of the hadron 
spectrum at high Pt does not come from thermalized matter but from QCD jets, yet it 
bears the imprint of the quark gluon plasma that the particles had to traverse in order to 
arrive at the detectors. It had already been realized in the 1980s that we can use high-Pr 
hadrons as "hard probes" of the quark gluon plasma phase [7]. 

In this overview we briefly cover some of the most fundamental approaches to heavy 
ion collisions from the theoretical side. We will discuss initial interactions and the color 
glass condensate, the current state of the lattice QCD results, jet quenching, and photons 
as an example for electromagnetic probes. We have to forgo some other very timely and 
important topics like heavy quarks, dileptons or hydrodynamics. The latter wiU be covered 
in an accompanying article by R. Bhalerao [8]. 



2. First Interactions 



Interactions of hadrons or nuclei at small momentum transfers, below 1-2 GeV/c can not 
be described by standard perturbative QCD (pQCD) techniques. pQCD assumes rather 
dilute parton densities in the colliding hadrons and is vahd at weak coupling. Experimen- 
tal data show that the gluon distribution grows dramatically for larger energies (or small 
parton momentum fraction x) which must eventually lead to a breakdown of the standard 
pQCD picture. At very large energies the distribution of gluons in hadrons (and even more 
so in large nuclei) has to saturate, with gluon splittings and recombination happening in 
detailed balance. It has first been argued by McLerran and Venugopalan that in such a 
saturation scenario the gluon density ~ Q sets a new scale, the saturation scale, and that 
an effective theory can be constructed based on a quasi-classical approximation to QCD 
[9-11]. For Qs ^ Aqcd the coupling can still be weak, although the theory has non- 
trivial many-body effects included. This effective theory of QCD at very high collision 
energies, called the color glass condensate, describes a nuclear collision as the colUsion of 
two (highly Lorentz-contracted) sheets of color S'C/(3)-charges [12,13]. 

The initial gluon fields Aj*, At^ in the two nuclei before the collision are dominated 
by their transverse electric and magnetic components. After the collision, this creates 
extremely strong longitudinal (color-)electric and magnetic fields [14,15] 

Eo = tg[A\,Al], B„ = ige''[A\,Ai] (1) 

(expressed here in a suitable axial gauge), see Fig. 1. Note that the nuclei themselves 
(represented by the conserved net baryon number carried by the valence quarks) go through 
each other Valence quarks (and other large-a; partons) only scattering occasionally to 
create jets. The longitudinal fields expand in the space between the nuclei as they recede 
from each other. We now understand that these strong longitudinal fields are the means by 
which energy density is deposited in the center of the colliding system. The initial energy 
density is given by 

- 2 l-^o + ^o) ~ in -J^ (2) 

where the /if are the densities of color charges in the two colUding nuclei and Qq and A 
are ultraviolet and infrared cutoffs. 
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Figure 1. Two Lorentz-contracted nuclei shortly after the collision. The large momen- 
tum partons have gone through each other with small disturbance. Their gluon fields 
have interacted to form strong longitudinal fields which connect the receding nuclei and 
slow them down. The scenario is reminiscent of a "color" capacitor with moving plates. 

The initial energy momentum tensor is diagonal with maximum pressure anisotropy. 
The "transverse pressure" — a slight abuse of terms so far from thermal equilibrium — is 
prp = T^^ = T^y — ep and the "longitudinal pressure" is pL = T^^ = —pT = — eo- The 
negative longitudinal pressure corresponds to the fact that the fields slow down the sources 
and hence the nuclei themselves (similar to the electric field attracting two capacitor plates 
in electrodynamics). The strong longitudinal fields lead to the creation of transverse fields 
as time progresses and the anisotropy between transverse and longitudinal components of 
the energy momentum tensor decreases. Equilibration needs two important ingredients: 
the gluon field has to decohere into gluons and quark-antiquark pairs in order to make the 
longitudinal pressure positive, and eventually the pressure has to isotropize completely. 
The exact way in which thermalization into an equilibrated quark gluon plasma proceeds, 
at least within the very short time (j 1 fm/c) suggested by data, is unknown [16,17]. Plasma 
instabilities have recently been suggested as a possible mechanism [18,19]. The time for 
decoherence, a necessary ingredient for thermalization, was recently estimated to be pro- 
portional to 1/Qs with a coefficient of order 1 [20,21]. For RHIC energies where Qs is 
estimated to be in the range 1-2 GeV, decoherence times are thus consistent with the exper- 
imentally observed thermalization times. Once thermalization is reached, the expansion 
and cooling of quark gluon plasma is very well described by hydrodynamics. The strong 
longitudinal gluon fields could play an important role in two recent discoveries at RHIC: 
long-range rapidity correlations which could be a direct image of the elongated flux tubes 
[22], and parity violation from quarks interacting with topologically non-trivial gluon con- 
figurations associated with E • B 7^ [23]. 

3. Lattice Estimates 

In lattice QCD the partition function of QCD is evaluated on a discretized Euclidean space- 
time grid using Monte-Carlo techniques. Results have long confirmed our expectation, 
originally from considerations of asymptotic freedom, that quarks and gluons are decon- 
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fined at high temperatures and that chiral symmetry is restored. There are strong indica- 
tions from reahstic simulations of QCD with two light and one heavier quark flavors (for u, 
d, and s quarks resp.) that the phase transition is not of first order at vanishing net baryon 
densities, but a smooth cross over around a critical temperature Tg. Hence we do not ex- 
pect a sharp transition, say of the energy density, but rather the thermodynamic properties 
change rapidly but smoothly as a function of temperature around Tc. 

Lattice QCD at finite baryon chemical potential /xb suffers from an unpleasant sign 
problem that makes Monte-Carlo techniques rather ineffective. In recent years, several 
techniques have been devised (reweighting, Taylor expansion around /ig = 0, etc.) to 
explore the phase diagram of QCD away from zero net baryon density. Most calculations 
agree that there is a critical point in a range /zb~200...400 MeV from which on the cross 
over line between hadronic and quark gluon matter becomes a first order phase transition. 
Our understanding of the QCD phase diagram is shown schematically in Fig. 2. 




Figure 2. Schematic phase diagram of QCD. The first-order phase transition Une ends 
in a critical point, at smaller baryon densities the transition between hadronic matter 
and quark gluon plasma is a cross over. 

Many modern calculations use close to realistic light quark masses with pion masses as 
small as 200 MeV. Estimates for the critical temperature at small baryon densities — the 
domain in which RHIC at top energies and the Large Hadron Collider (LHC) are operating 
— range from 176 to 192 MeV [24,25]. This is a rather large discrepancy for results 
from different groups. While many results hint at one common temperature for the chiral 
and deconfinement phase transition, one group sees evidence for a chiral transition at 151 
MeV while they measure the deconfinement transition at 176 MeV [24]. The divergence 
of predictions has recently attracted a lot of attention and is under investigation. 

There is however consensus on the gross features of the QCD equation of state calculated 
on the lattice. Around the critical temperature the normalized energy density e/T^ exhibits 
a large jump due to the drastic increase in degrees of freedom going from the hadronic 
phase (basically pions below Tc) to the deconfined phase with quarks and gluons. e/T^ and 
the reduced pressure p/T^ do not reach their Stefan-Boltzmann limit even around 4 x 
which has been interpreted as a hint for the strongly interacting nature of quark gluon 
plasma even far above the phase transition [4]. However, refined perturbative calculations 
can describe the lattice results reasonably well down to about 2 x Tc [26]. 

More recently it has been attempted to extract transport coefficients like the shear vis- 
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cosity 1] and the bulk viscosity of quark gluon plasma from lattice QCD calculations. As 
improved viscous relativistic hydrodynamic calculations become more and more available 
the interest in these new observables will increase [8]. 



4. Hard Probes 



Jet quenching, the loss of high momentum particles in nuclear collisions had long been 
predicted theoretically [7]. When it was finally observed in early RHIC data the effect 
turned out to be very large. Nearly 80% of the particles expected at a given transverse 
momentum Pt were suppressed in central collisions of gold nuclei. We quantify this 
suppression most conveniently by the nuclear modification factor 

{N,on)dNPP/dPT ^ ' 

which is the ratio of particle yields in nucleus-nucleus vs proton-proton collisions modulo 
the number of binary collisions expected in a nuclear collision. For a loose superposition of 
nucleon collision the modification factor would be (close to) unity. In fact there are a few 
"cold" nuclear matter effects which lead to deviations of Raa from unity. These include 
modifications of parton distributions for nucleons bound in nuclei vs free nucleons (known 
as shadowing and the EMC effect), and multiple scattering in the initial state which leads to 
the so-called Cronin effect [27]. At top RHIC energies we expect particles produced with 
a few GeV energy around midrapidity to exhibit anti-shadowing, i.e. a slight enhancement, 
which gives way to a suppression at larger Pt- The Cronin effect leads to an enhancement 
for the same particles at a few GeV which dies out toward higher Pt. However, all of these 
cold nuclear matter effects are small, and partially competing with each other. In summary 
they lead to deviations of Raa from unity of less than 20%. 

Two effects are important to understand the energy loss of a fast quark or gluon prop- 
agating through quark gluon plasma. First, the dominating mechanism for energy loss, at 
least for light quarks and gluons, is through induced gluon radiation. This leads to larger 
suppression than energy loss than from elastic collisions. Secondly, the dependence of the 
energy loss on the path length L is quadratic, AE oc L^. This LPM effect (after Landau, 
Pomeranchuk and Migdal who first described the similar effect in quantum electrodynam- 
ics [28,29]) comes from a destructive interference. Consider a parton radiating a gluon 
with energy lj and relative transverse momentum kx- The formation time of the parton- 
gluon system is given by r/ = ui/k^, and further radiation during this time is suppressed. 
In other words, if the mean free path A of the parton is smaller than the typical formation 
time Lj/k^, then the scattering off A^coh = ^z/A number of partons happens coherently. 
This leads to a differential energy loss dE/dx = —qx per unit path length, and to the 
famous quadratic path dependence of the total energy loss. The transport coefficient q pa- 
rameterizes the momentum transfer between the medium and the parton. It is given by the 
momentum transfer squared /x^ per mean free path. 

Comparisons of q extracted from measurements at RHIC and from cold nuclear matter 
experiments, e.g. HERMES, show a large increase of the quenching power at RHIC [30]. 
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The handwaving arguments above can be backed up by calculations which compute 
energy loss under certain simplifying assumptions. The most important models based on 
perturbative QCD are 

• the Higher Twist (HT) formalism by Guo and Wang [31,32,30]. It uses results from a 

rigorous calculation of medium-modified fragmentation functions in semi-inclusive 
deep-inelastic electron scattering off nuclei (e+A) and transfers these fragmentation 
functions to nuclear collisions. 

• the AMY formalism based on a series of papers by Arnold, Moore and Yaffe [33,34]. 
These authors use hard thermal loop resummed perturbation theory, valid at asymp- 
totically large temperatures, to calculate the rate of energy loss for fast partons 

• the ASW formalism developed by Armesto, Salgado and Wiedemann [35,36]. Here 
multiple soft gluon emission is resummed using a Poisson statistics for the number 
of gluons emitted. 

• The GLV approach by Gyulassy, Levai and Vitev [37-39]. They discuss scatterings 
off static scattering centers in an opacity expansion in the medium. 



As an example we briefly highlight some fundamental ideas of the Higher Twist for- 
malism. Fig. 3 shows typical diagrams contributing to the usual DGLAP equations for 
fragmentation functions (left panel), and to a set of modified evolution equations that are 
formally of higher twist in deep-inelastic scattering (DIS). Higher twist in factorized per- 
turbation theory usually refers to sub-leading contributions that are suppressed by powers 
of the ratio of the non-perturbative over the large perturbative scale, ^ A^/Q^ << 1. In 
nuclei, some of those higher twist contributions are re-enhanced through a size factor, such 
that they enter at a level LA^/Q"^ ~ 1. Note that we have indicated two propagators each 
in amplitude and complex conjugated amplitude in the right panel through grey and black 
circles. The interference between the poles of two of those propagators leads to the LPM 
suppression effect in the Higher Twist formalism. The modified evolution equations define 
medium-modified fragmentation functions D{z) which describe the energy loss of partons 
in semi-exclusive DIS. In their simplest form they can be parameterized as 




Figure 3. Left panel: diagram contributing to the DGLAP evolution equation of (vac- 
uum) fragmentation functions in semi-exclusive deep-inelastic scattering. Right panel: 
typical diagram contributing to the modified evolution equations at the twist-4 level. 
The outgoing parton scatters off the medium and radiates. 




(5) 
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where z is the momentum fraction of the hadron within the fragmenting parton and Az 
is the average shift due to induced radiation. Similar fragmentation functions, with the 
average energy loss scaled up to fit the data, can then be used to describe high-Pr data at 
RHIC. 

The energy loss models described above differ in some of their underlying assumptions 
and are by no means completely equivalent. Nevertheless, all of them manage to describe 
both the Py-dependence and the impact parameter-dependence of Raa measured at RHIC 
remarkably well after one parameter, typically ^ or a related quantity, is fixed. With the 
advent of data on hadron correlations at high Pt additional non-trivial restrictions on fits 
were available. Typically experiments measure pairs of hadrons, a trigger hadron t and an 
associated hadron a which are then presented as a per-trigger yield 

dN/dPtdP^d{A<j>) 

F(P„P.,A0) = (6) 

as a function of the transverse momenta of the trigger and associate particle, Pt and Pa 
resp., and the relative azimuthal angle Acj). To determine the suppression relative to free 
proton-proton collisions we can once more look at the appropriate nuclear modification 
factor 

^ Y^^{PuPa,A^) 

AA YPP{Pt,Pa,A(P) ■ 

Iaa poses the first true challenge for jet quenching models, and some studies have found 
incompatible values of q from separate fits to Raa and Iaa- It has also become clear that 
there is a long list of uncertainties, e.g. details of the modeling of the background fireball 
[40], the treatment of energy loss at early times before the formation of equilibrated quark 
gluon plasma, etc. [41]. 

Where do we stand with hard probes? A comparative study by Bass et al. [42] found 
unacceptably large discrepancies between values of q derived from different jet quenching 
models fitted to RHIC data. The extreme values found in this study are 

g = 18.5 GeVVfm for ASW , q = 4.5 GeV Vfm for HT . (8) 

The good news is that these values are, by far, larger than for cold nuclear matter. How- 
ever, they are not compatible among each other. A sustained effort is now underway, 
spear-headed by the TECHQM initiative and the JET collaboration, to systematically in- 
vestigate hard probes, by comparing and vetting different calculations, and by identifying 
and narrowing down the sources of theoretical uncertainties. 

Once a reliable value of q is established, a comparison with its perturbative value in an 
equilibrated plasma Qpcrt = 2e^/^ will tell us whether the coupling of jets of the quark 
gluon plasma is weak or strong. Some other interesting developments are the emergence 
of simulations for the full evolution of jet showers in a medium [43,44], medium modifi- 
cations to jet shapes [45], and hadron chemistry at large momentum Pt [46]. 



5. Photons 



Hadrons from the bulk carry information from the point of their last interaction (the freeze- 
out), and hadrons at large momentum Pt carry some information of the fireball integrated 
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over its history. On the other hand electromagnetic probes, i.e. real and virtual photons 
(dileptons) have the wonderful property that they do not re-interact once created, due to 
their large mean-free path in quark gluon plasma. Hence they are unique penetrating probes 
that can give us unobstructed information about the center of the fireball in nuclear colU- 
sions, and about the earUest times during the collision. The following sources of photons 
are important in relativistic heavy ion collisions [47]. 

• Prompt hard photons from hard collisions: they are created from annihilation and 
Compton scatterings, q + q ^ j + g and g + 5 — > 7 + 3 resp., of large momentum 
partons in the original nuclear wave functions. They dominate the direct photon 
spectrum at the largest transverse momenta Pt- They also do not carry information 
about the quark gluon plasma formed, but exhibit some of the initial state effects 
(shadowing, Cronin effect) mentioned above. 

• Fragmentation photons from vacuum bremsstrahlung: any initial hard scattering of 
partons can lead to bremsstrahlung. Most of the time this will be gluon radiation, 
but there is a sizable contribution of photons which can be described by parton-to- 
photon fragmentation functions. Vacuum fragmentation photons have a somewhat 
steeper _Pt -spectrum but are comparable in yield with prompt hard photons except 
for the largest momenta. Their yield can be modified in nuclear collisions as partons 
travel over some distance and lose energy before radiating a photon. 

• Photons from jet conversions and induced bremsstrahlung [48-51]: as leading jet 
partons travel through the medium and lose energy through induced gluon radiation, 
again part of this induced radiation can be in the form of photons. In addition, elastic 
annihilation and Compton scatterings of fast jet partons with partons in the quark 
gluon plasma can lead to an effective conversion in which all of the momentum of 
the fast parton can be transferred to the photon. These processes give competitive 
photon yields in heavy ion collisions at RHIC and LHC at intermediate Pt of a few 
GeV/c. The photons carry the full time evolution of the quark gluon plasma fireball 
as the local conversion rates are proportional to In T and thus one can estimate 
plasma temperatures T from their yields. 

• Thermal photons from both the QGP and the hadronic phase [52,53]: annihilation 
and Compton as well as bremsstrahlung processes between quarks and gluons also 
play out in the equilibrated plasma, providing an exponential photon spectrum that 
dominates at low Pt- Its slope obviously provides the temperature T, although the 
photons seen in the detector are integrated over the spatial profile as well as the 
time evolution of the fireball. The hot hadronic gas phase below the phase transition 
temperature also leads to thermal photons at lower temperatures. 

All of these sources have been computed in a sustained effort over two decades. A mod- 
ern, comprehensive calculation and comparison with data can be found in publications of 
the McGill group [54]. They compare well with the experimentally measured spectra and 
Raa of direct photons (decay photons of hadrons like the tt" bury the direct photon signal 
in experimental data and have to be subtracted). At low Pt there is ample evidence that the 
"glow" of the quark gluon plasma has been found and that the temperatures are above the 
predicted values of Tc- Constraints on other sources, like the conversion photons, are not 
as stringent. Other observables have been suggested to explicitly find those photons, and 
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to get a complementary picture of both the medium and jet quenching: e.g. the azimuthal 
asymmetry t'2 of photons, and photon-hadron correlations. 

V2 is the second order harmonic of the spectrum dN/ d(j) w.r.t. the azimuthal angle around 
the beam axis. The reference axis (0 = 0) is given by the reaction plane. So far we have not 
discussed V2 for hadrons from the bulk fireball or for hadrons from jets. In both cases V2 is 
positive, meaning more particles are emitted into the reaction plane than out of it. For bulk 
particles this comes from the hydrodynamic expansion with larger pressure gradients and 
hence larger flow into the plane, for high- Ft particles it derives from a smaller opacity and 
less absorption into the plane. It has been predicted in [55] that more conversion photons 
should be emitted in the elongated direction of the fireball, out of the reaction plane than 
into the reaction plane. This would lead to a negative V2 for this photon source, the first 
process discovered to exhibit this behavior Its experimental discovery would unanimously 
confirm conversion and induced bremsstrahlung of jets in medium. Since other photon 
sources come with vanishing or positive V2 the sign of the total direct photon V2 is rather 
unclear and the absolute value is small due to the cancellations. Experimental data confirm 
small absolute values but the error bars in the Pt region of interest are too small to pin 
down the sign [55,56|. 

One of the most exciting tools in our arsenal combine properties of hard and electromag- 
netic probes.These are photon-hadron correlations at high momentum. Photon-triggered 
fragmentation functions can give unprecedented access to the mechanisms of energy loss, 
as prompt hard photons carry the information about the initial momentum of the hard recoil 
quark or gluon. This information is not available if Raa is considered and it is consider- 
ably washed out even if Iaa is considered. First experimental results of photon-triggered 
fragmentation functions are available. 

6. Summary 

The processes unfolding in collisions of nuclei at high energies are qualitatively under- 
stood. The system goes from a phase of strong initial gluon fields through a thermal- 
ization process to a hydrodynamically expanding quark gluon plasma which eventually 
hadronizes. QCD jets, high-Px hadrons and photons can serve as built-in probes of the 
new phases of QCD. The field is now moving into a stage of quantitative analysis in which 
transport coefficients and other observables should be extracted with much reduced error 
bars. In the near future the high energy ion program at LHC and the future FAIR facility 
wiU provide completely new perspectives on the QCD phase diagram. 
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